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We present a detailed description of how a differential geometric approach to 
Hamiltonian dynamics can be used for determining the existence of a crossover 
between different dynamical regimes in a realistic system, a model of a rare 
gas solid. Such a geometric approach allows to locate the energy threshold 
between weakly and strongly chaotic regimes, and to estimate the largest 
Lyapunov exponent. We show how standard methods of classical statistical 
mechanics, i.e. Monte Carlo simulations, can be used for our computational 
purposes. Finally we consider a Lennard Jones crystal modeling solid Xenon. 
The value of the energy threshold turns out to be in excellent agreement 
with the numerical estimate based on the crossover between slow and fast 
relaxation to equilibrium obtained in a previous work by molecular dynamics 
simulations. 
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I. INTRODUCTION 



Generic non-integrable Hamiltonian systems with N > 3 degrees of freedom always have 
a connected chaotic component in phase space. Moreover, as N is large, the measure of 
such a component should be practically coincident with the measure of the whole constant- 
energy hyper-surface. In fact, as discussed in Refs. ||T]-§|, the invariant tori whose existence is 
predicted by the Kol'mogorov-Arnol'd-Moser (KAM) theorem have a positive measure only 
below a critical value for the amplitude of the non-integrable part of the Hamiltonian, and 
this critical amplitude is estimated to be rapidly decreasing with N. This does not mean 
that particular cases in which the KAM threshold is relevant also in macroscopic systems 
cannot exist, nevertheless it suggests that such a situation could be hardly generic. 

These facts support the expectation that varying the energy - or more precisely the 
energy density e = E/N which is the physical parameter as N is large and eventually as the 
limit N — > oo is taken - one should not observe any qualitative change in the dynamical 
behavior of large non-integrable Hamiltonian systems. The dynamics should be completely 
chaotic at all energies, and the only effect of a variation of the energy should be a somewhat 
trivial rescaling of the characteristic instability time scale, measured by the inverse of the 
Lyapunov exponent A. An example of this kind of behavior is provided by self-gravitating 
systems [|J. 

On the contrary, there is now a widely accepted numerical evidence that - at least as long 
as models of nonlinear coupled oscillators are considered - there exist qualitatively different 
regimes in the dynamics, which have been referred to as weak and strong chaos Jl],|2| . In the 
strongly chaotic regime which corresponds to the above sketched scenario, fast phase-space 
mixing is observed regardless of the initial conditions. At variance, in correspondence of 
weak chaos one can observe very long mixing times with non-equilibrium initial conditions, 
and the details of the dynamics are strongly influenced by the choice of the initial conditions. 
Moreover, at least on finite time scales, the dynamics appears as globally recurrent, as is 
probed by the probability distribution of single-particle autocorrelation functions ||. This 
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effect was observed in several numerical simulations |||7]] inspired by the results of the 
celebrated numerical experiment by Fermi, Pasta and Ulam (FPU) || where the expected 
equipartition of energy among normal modes was not observed in a chain of linear oscillators 
coupled by a weak anharmonicity. 

The transition between weak and strong chaos is rather sharp when detected looking at 
non-equilibrium properties, i.e. observing the time behavior of observables which depend 
on the choice of a particular dynamical initial condition in which the system is far from 
thermodynamic equilibrium. Examples of these observables are the relaxation times of 
several dynamical observables to their equilibrium value |l|J^] , the finite-time values of the 
so-called spectral entropy |7]|§, and the already mentioned probability distributions of single- 
particle autocorrelation functions ||. The behavior of these observables neatly detects a 
threshold value e = e c which marks the transition between weak and strong chaos, and has 
been referred to as the strong stochasticity threshold (SST) fl|-f|,|T0|| , or as the crossover 



energy (CE) |TI[. All these observables have the drawback of being not globally defined, i.e. 
of depending on a particular choice of the initial condition, which could in principle depend 
on N. 

It is remarkable that also the Lyapunov exponent A, an observable which is neither 
dependent on the initial conditions nor on N, and which measures directly the degree of 
chaos, marks the SST. In fact in correspondence of e c the dependence of A on e has a 
crossover: in the strongly chaotic regime one finds a power law which can be successfully 
predicted by a random matrix approximation (RMA) for the tangent dynamics |J. In 
this regime the RMA yields good predictions also for the shape of the whole spectrum of 
Lyapunov exponents fl2| . At variance, in the weakly chaotic case A is still positive - for 



nonlinear coupled oscillators A oc e 2 - thus indicating that chaos is present, but the RMA 
is no longer able to predict its dependence on e. The breakdown of the RMA is a clue 
that some global change has happened in the phase space structure. The drawback of using 
the dependence of A on the energy density as a probe of the SST is that the transition is 
no longer sharp as in the case of non-equilibrium observables: at e c , A(e) exhibits only a 



crossover between different asymptotic behaviors. 

The deep origin of these puzzling dynamical features of non-integrable Hamiltonian sys- 
tems with many degrees of freedom has at present not yet been understood: nevertheless, a 
recently proposed differential geometrical approach to Hamiltonian chaos H has established 
a link between the SST and some major change of the geometric structure underlying the 
dynamics, allowing an operative definition of this threshold no longer based on the compu- 
tation of time-asymptotic quantities like Lyapunov exponents, but on statistical averages of 
geometric observables. Starting from the results of Ref. || it has been proved the stability of 
the SST in the thermodynamic limit N — ► oo, at least for the one- dimensional FPU model 
T0|| • Moreover this approach provides the basic tools to obtain a model scalar equation 



which describes the main features of chaos being independent of the details of the dynamics 
|T3| . The present work follows this geometric approach. 

It is worth mentioning that very recently another kind of geometric approach has been 



put forward jT4| which is based on the geometric properties of the trajectories - seen as 
curves in M? N - rather than on the properties of the ambient manifolds. As long as a 
comparison is possible, the two approaches yield perfectly consistent results. 

The phenomenology related to the SST has been mainly studied in connection with 
one-dimensional models of coupled oscillators (FPU model and lattice if 4 classical model 
JTJ-^J]) , and some numerical results for two-dimensional crystals [jljj strongly suggest that a 
transition between weak and strong chaos is present in more realistic systems. The present 
work is concerned with the possibility of the existence of the SST in a realistic system, in 
the perspective of a possible experimental verification of the physical consequences of the 
transition between weak and strong chaos. Such an experiment has been recently proposed 



TT|j , and should verify the existence of a crossover between slow and fast relaxation to 
thermal equilibrium in a rare gas crystal with diluted impurities. The far-from-equilibrium 
dynamics of such a system has been studied numerically in Ref. |TTJ, showing that the 
crossover occurs at e ~ 0.15s, where e is the depth of the interaction potential well. The 
present work is concerned with the equilibrium dynamical properties of a simpler but closely 



related system, i.e. a model of a rare gas solid without impurities. As discussed in Ref. |TT 



the role played by impurities is crucial to obtain suitable initial conditions in the non- 
equilibrium case, but should have only a weak effect on the global dynamical properties of 
the system (this fact is confirmed in a 1-d case |16|]), hence the eventual existence of a SST 
in our model system - apart from being interesting by itself - could have a significance in 
explaining the phenomenology observed in Ref. ]ll] . 

More precisely, the model studied in the present work is a system of n point masses 
arranged on a three-dimensional fee lattice with nearest-neighbor Lennard- Jones interaction. 
The classical dynamics of such a system is a reasonable approximation of the behavior of 
a rare gas crystal if quantum effects can be neglected, i.e. as long as the temperature is 
high enough. Hence the existence of a crossover between weak and strong chaos in this 
model could have a substantial physical relevance if it occurs at e values allowing a classical 



description of the dynamics, at least as a first approximation. As we shall see in Sec. |TV], 
this happens at least in the case of Xe crystals. 

The present work is organized as follows: in Sec. [TT] the geometrical methods which 
allow a characterization of the transition between weak and strong chaos are sketched; these 
methods were introduced in Refs. P,|TU| and where all details can be found. The model 



studied is described in Sec. Ill], and the results are presented and discussed in Sec. |TV 
Some conclusions are drawn in Sec. [V|. 



II. RIEMANNIAN GEOMETRY AND CHAOTIC DYNAMICS 

Hamiltonian dynamics can be rephrased in geometrical terms owing to the fact that the 
trajectories of a dynamical system with quadratic kinetic energy can be seen as geodesies 
of a suitable Riemannian manifold. There are several choices for the ambient manifold 
as well as for the metric tensor. As already discussed in Ref. QTUJ a particularly useful 
ambient space is the enlarged configuration space-time Mxl 2 , i.e. the configuration space 
{g 1 , . . . , q\ . . . , q N } with two additional real coordinates q° and q N+1 . In the following q° will 
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be identified with the time t. For standard Hamiltonians Tt = T + V(q) where T = ha>ijQ l $, 
this manifold, equipped with Eisenhart's metric gE, has a semi-Riemannian (Lorentzian) 
structure (detg^ = — 1). The arc-length is given by 

ds 2 = UijdqW - 2V(q){dq ) 2 + 2dq°dq N+1 , (1) 

where both i and j run between 1 and N. Let us restrict to geodesies whose arc-length 
parametrization is affine, i.e. ds 2 = 2C 2 dt 2 ; simple algebra shows that the geodesic equations 
d=Y , r/1 dq v dq x 

become Newton equations (without loss of generality = 6ij is considered) 

d 2 q { _dV 
~d¥ ~ ~dq~ 

for % = 1 . . . N, together with two extra equations for q° and q N+1 which can be integrated 
to yield 

q° = t (4a) 



q N +i = C\t + C 2 - I L(q, q) dt (4b) 



o 

where L(q, q) is the Lagrangian, and C\, Gi are real constants. As stated by Eisenhart 
theorem the dynamical trajectories in configuration space are projections on M of the 
geodesies of (M x R 2 ,g E ). 

In the geometrical framework, the stability of the trajectories is mapped on the stability 
of the geodesies, hence it can be studied by the Jacobi equation for geodesic deviation 

^ + fl(7, J) 7 = , (5) 

where R is the Riemann curvature tensor, 7 is the velocity vector along the reference geodesic 
7(5), D/ds is the covariant derivative and J, which measures the deviation between nearby 
geodesies, is referred to as the Jacobi field. The stability - or instability - of the dynam- 
ics, and thus deterministic chaos, originates from the curvature properties of the ambient 
manifold. In local coordinates, Eq. |5| is written as 



and as already shown in Refs. [|3|,|T0ll , in the case of Eisenhart metric it simplifies to 



d 2 J l d 2 V 

~w + d^ J3 = ' (7) 

which is nothing but the usual tangent dynamics equation for standard Hamiltonians. The 
Lyapunov exponents are usually computed evaluating the rate of exponential growth of J 
by means of a numerical integration of Eq. ([?D |Tj| . 



In the particular case of constant curvature manifolds, Eq. @ becomes very simple JT9 

D 2 

and has bounded oscillating solutions J ~ cos(V^s) or exponentially unstable solutions 



J pa exp(\/—K s) according to the sign of the constant sectional curvature K, which is given 
by 

K = J^= U (9) 
N-l N(N — 1) ' [ ) 

where K R = R fiu q fi q u is the Ricci curvature and 71 = is the scalar curvature; R^ u is the 
Ricci tensor. Manifolds with K < are considered in abstract ergodic theory (see e.g. Ref. 
). Krylov [21] originally proposed that the presence of some negative curvature could 



be the mechanism actually at work to make chaos in physical systems, but in realistic cases 
the curvatures are neither found constant nor everywhere negative, and the straightforward 
approach based on Eq. @ does not apply. This is the main reason why Krylov' s ideas 
remained confined to abstract ergodic theory with few exceptions. 

In spite of these major problems, some approximations on Eq. (|5|) are possible even in 
the general case. The key point is that negative curvatures are not strictly necessary to 
make chaos, and that a subtler mechanism related to the bumpiness of the ambient manifold 
is actually at work. 

Let us choose a geodesic frame (i.e. a reference frame which is parallel transported along 
a geodesic; as a consequence, D/ds = d/ds), and project Eq. (|5|) on a direction determined 
by the unit vector u: 



^(J,u) + (R(i,J)i,u) = 0. (10) 

If the system is chaotic, J grows exponentially with growth rates given by the Lyapunov 
exponents, and if u is the direction corresponding to the largest Lyapunov exponent A, 
after a finite (proper) time s the components of J along the other directions will become 
negligible compared to that along u, thus we find J ~ ipu. Equation flip]) is thus rewritten 
approximately as a scalar Hill equation for if), 

^ + K(s)ij = 0, (11) 

where K(s) = K(j, u) is the sectional curvature of the geodesic plane spanned by the direc- 
tions 7 and u and is no longer a constant, but a fluctuating function taking mostly positive 
values (in some cases like the FPU model K is strictly positive) whence the solutions of Eq. 
(|TT|) can be subject to parametric instability. Curvature fluctuations can produce chaos even 
if no negative curvature is experienced by the geodesies. As the sectional curvature is no 
longer constant, Kr and 1Z are respectively averages of K over the direction of J and over 
both the direction of J and the direction of the reference geodesic in the latter. Equation 
(§) no longer holds, nevertheless it is a first order approximation to which an estimate of the 
curvature fluctuations can be added to obtain a stochastic model of K(s) independent of the 
dynamics of the system . This model leads to an analytical estimate of the Lyapunov 
exponent A, which is correct (at least for the FPU model) in the limit iV — > oo. 

Up to this point the results are independent of the choice of the metric. Specializing to 
the Eisenhart arc- length parametrization Eq. (|11|) is rewritten in terms of the time t, 

i> + K{t)ip = , (12) 

where a dot stands for a time derivative and K{t) = K(s^/2Cf) /2Cf. The stochastic model 
of K(t) is given by 

K(t) = (k R ) + (Pkn) 1 ' 3 *® , (13) 
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where kn = Kr/N, (•) stands for an average taken along a geodesic, which, for systems in 
thermal equilibrium, can be substituted with a statistical average taken with respect to a 
suitable probability measure (e.g. the micro-canonical or the canonical measure); rj(t) is a 
stationary 5-correlated Gaussian stochastic process with zero mean and variance equal to 
one. Using Eisenhart metric, and for standard Hamiltonians, the non-vanishing components 
of the Riemann tensor are RoiOj = d qi d q .V, hence the Ricci curvature has the remarkably 
simple form 

k R = i W , (14) 
where V 2 is the Euclidean Laplacian operator. Equation ( |T^ ) becomes a stochastic differ- 



ential equation, i.e. the evolution equation of a random oscillator ||22| . It is worth noticing 
that Eq. (|1^) is no longer dependent on the dynamics, since the random process depends 
only on statistical averages. The estimate of the Lyapunov exponent A is then obtained 
through the evolution of the second moments of the solution of ( |l2|) as 

A . lim i log i!w±fw . (15) 

^oo2 6 ^ 2 (0)+V 2 (0) 



1 / 4fc\ 

\(k,5 k ,T) = -[A- — ) , (16) 



As shown in Ref. fL3| , this yields the following expression for A: 

2 V 3a; 

where 

1/3 



A = K^+\/^ + ^ 2 i . (17a) 
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7T 



(17b) 



2y/k(k + 6 k )+Tr6 k 

in the above expressions k is the average Ricci curvature k = (kji) and 5 k stands for the 
mean-square fluctuation of the Ricci curvature, 5 k = (5 2 kR) 1 ^ 2 . 

The advantages in using the geometric approach to Hamiltonian chaos are thus evident. 
In fact, it is possible to give reliable estimates of the Lyapunov exponent without actually 



computing the time evolution of the system: the estimate (|16D of A depends only on statistical 
averages which can be either computed analytically in some cases (for instance in the case 
of the FPU model or, in general, extracted from a Monte Carlo simulation, as it is the 
case of the model to be studied in the present work. 

The behavior of the average geometric observables as the control parameter (e.g. the 
energy density or the temperature) is varied conveys an information which goes beyond 
the possibility of computing the Lyapunov exponent. The dependence of the average Ricci 
curvature on the energy density has already been used in Ref. [Kj to give an operational 



definition of the SST which allows its computation in the thermodynamic limit, showing the 
stability of the threshold in this limit for the FPU model. In fact, it is easy to find that 
for a harmonic chain (k^) is constant as the energy density is varied, and this is a common 
feature of other integrable models (e. g. the Toda chain |10[)- A computation of {k R )(e) at 
constant volume (length) for the FPU chain shows that the average Ricci curvature exhibits 
two well-defined asymptotic behaviors, {k R ){e) = const, as e — > (harmonic limit), and 
(kii)(e) = e 1//2 as e — > oo. The crossover between the two asymptotic curves occurs at a 
value of e, e c , which can be interpreted as a geometric estimate of the SST. Such a geometric 
estimate is in very good agreement with estimates based on other methods. 

Moreover, one can look at the random oscillator equation fll2l ) as an effective Jacobi 
equation for a geodesic flow on a surface M whose Gaussian curvature is given by the 
random process K(t). As long as nonlinear coupled oscillators are considered, the average 
Ricci curvature is positive, hence M can be regarded as a sphere with a fluctuating radius. 
In the limit of vanishing fluctuations, one recovers the bounded evolution of the Jacobi field 
associated with integrable dynamics. Chaos suddenly appears as curvature fluctuations 
are turned on, nevertheless it it will be "weak" as long as 5k <C k, i.e. as long as M 
can be considered as a weakly perturbed sphere. On the contrary as the size of curvature 
fluctuations becomes of the same order of the average curvature, 5k ~ k, M can no longer 
resemble a sphere, and the dynamics will no longer "feel" the integrable limit. Hence we 
expect the dynamics to be strongly chaotic. This is by no means a deep explanation of the 
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existence of weakly and strongly chaotic regimes in Hamiltonian dynamics. Nevertheless 
it shows how the simple geometric concepts which enter the Riemannian description of 
Hamiltonian chaos, besides providing effective computational tools, are also useful in helping 
one's physical intuition with images and analogies which would be difficult to find elsewhere. 

III. THE MODEL AND THE GEOMETRIC OBSERVABLES 

The system studied in the present work is a crystal of n atoms of mass m moving in 
three dimensions and interacting through a pairwise central potential v(r). Its Hamiltonian 
is 

n 

n = ^J2tf + V W (18) 

i=l 

where 

1 71 

V(X.) = -J2 v (\*i-*i\)i (19) 

i,i=l 

The geometric observables which, within the approximations described in the previous Sec- 
tion, are relevant to the dynamical instability of the system, are the average Ricci curvature 
of (M x M 2 , gs) and its fluctuations. In the case of a pairwise interaction potential the Ricci 
curvature turns out to be (see Eq. [L4|): 

i=i ii=i 11 JJ 

where it should be noticed that N = 3ra as it represents the number of degrees of freedom. 
The quantities which have to be determined are 

k = (k R ) (21a) 

Si = (S 2 k R ) = I ((4) - (k R f) . (21b) 

The probability measure which is usually employed for computational purposes is the canon- 
ical distribution so that the statistical averages of Eqs. can be written as 
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(f)((3) = Z c 1 J dX/(X)e-^ x ) , (22) 
where Zc is the configurational part of the canonical partition function, 

"dXe-^W. (23) 



As a matter of fact in the canonical statistical ensemble the role of control parameter is played 
by the inverse temperature (3 = 1/kgT where ks is the Boltzmann's constant. However, in 
the thermodynamic limit (N — > oo), the micro-canonical averages can also be obtained from 
the canonical ones. As long as the mean value k = (kit) is concerned, the canonical and 
microcanonical averages differ only by a 0(1/ N) correction. Hence in the thermodynamic 
limit 

(*s>W = (*R> , (24a) 

^=h-^h [lotZcm - (24b) 

As regards Sf. = (5 2 kn), one must keep in mind that fluctuations depend on the statistical 
ensemble. In fact the difference between canonical and microcanonical fluctuations does 
not vanish in the thermodynamic limit and the relation between these two quantities is, 
according to Ref. 

(5 2 k R ) (e) = (6 2 k R ) (0 (e)) + F ((3 (e)) , (25a) 



f^-^ffl , (25b, 



dp 

where e ((3) is given again by Eq . ( p4b|) and c v is the specific heat at constant volume. Thus 



the micro-canonical fluctuations can be obtained from the canonical ones provided that the 
values of the specific heat are known. 
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IV. RESULTS AND DISCUSSION 



The geometric observables described in the previous Sections have been evaluated for a 
Lennard- Jones face-centered cubic (fee) crystal whose Hamiltonian is given by Eqs. (|T8|) and 
(|19|) with a pair interaction potential which reads 



t>(r) = As 



a\ 12 (o 



r 



(26) 



Through an appropriate choice of the free parameters m (mass), e and a, this simple model 
is able to take account of most thermodynamical properties of rare gas solids. 

The statistical averages of Eqs. (|2lD have been caculated by means of a standard canoni- 
cal Monte Carlo algorithm where a simulation box of n = 256 particles subjected to periodic 
boundary conditions has been used. In all simulations only nearest neighbors interactions 
have been dynamically taken into account; the contributions of the interactions beyond the 
nearest neighbors shell has been considered in a "static approximation" in which the istan- 
taneous relative positions of the atoms are replaced by their equilibrium values. Apparently 
this procedure does not affect the evaluation of the geometrical observables; its advantage 
resides in the fact that allows one to employ the all neighbors parameters e and a which 
give a reasonnable representation of the real pair potential and the equation of state of rare 
gas solids p5| . 

All data reported here are given in dimensionless form by reducing them with respect to 
the "natural units" of the model; namely the reduced energy is measured in unit of e and 
the Ricci curvature in unit of e/a 2 as well as its mean-square fluctuation. 

We performed two distinct series of simulations. In the first series the density p of the 
crystal has been kept constant in order to compare the qualitative behaviour of the geometric 
observables with the known results for one- dimensional systems. In fact, computations of 
chains of anharmonic oscillators have been performed at constant length P, |10| , |T6"|| . The 



numerical outcomes of this series of simulations are reported in Figs. 1-3 where k, 5k/ k 
and the estimate A of the Lyapunov exponent computed according to Eq. (|I6"D are plotted 
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versus the energy per degree of freedom of the system. The additional data appearing in 
these figures refer to computations where the interaction potential (26) has been expanded 
up to the fourth order in the Taylor series. This procedure allows us to perform simulations 
at higher energies in order to provide a better representation of the crossover energy region. 
The results reported in Fig. 1 show that the SST, or crossover energy, can be located around 
e c ~ 0.15. This estimate of e c is confirmed by the results reported in Fig. 2, for Fig. 2 shows 
that in correspondence of this threshold the ratio 5k /k deviates remarkably from the low- 
energy behaviour 5k/ k oc e and tends to saturate towards 5k /k ~ 1 which implies that the 
manifold becomes highly anisotropic, and thus we expect the dynamics to enter the strongly 
chaotic regime. As a consequence of the behaviour of k and 5k, the geometric estimate A 
of the Lyapunov exponent shows a sharp crossover between two different power laws. It is 
worth recalling that the quantity here reported is not the "true" Lyapunov exponent but an 



estimate which has anyhow proven to be extremely accurate in other systems |T3 . 

The aim of the second series of simulations is mainly to determine whether the aver- 
ages of the geometric observables are affected by a change in the density of the sample. 
In principle one expects such a change because the curvature properties depend on space 
derivatives, and a change of the density induces a change in the length scale. Moreover, it 
is interesting to test the expectation that the crossover in the relaxation behaviour of a Xe 
crystal recently observed in numerical simulations ||11|| , is related with some major change 
in the global properties of the dynamics (i.e. the SST) which in turn can be detected by the 
geometric observables under investigation in this paper. Hence it is interesting to compute 
the statistical averages of Eqs. (|21| ) using settings which are close to those used in the 
molecular dynamics simulations reported in fLl| . In particular the results shown in Figs. 



3-6 refer to the values of the density given by the empirical equation of state of solid Xe 
(e/k B = 228.6° K, a = 3.959 A) @. 

The results reported in Figs. 3-6 show that the qualitative behaviour of the curvature 
fluctuations (Fig. 5) and of the theoretical estimate of the Lyapunov exponent A (Fig. 6) is 
the same as in the constant density case, while the behaviour of the average curvature (Fig. 
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4) is dramatically altered: instead of a crossover between two asymptotic regimes we have 
here a neat maximum of k at e = e c . It is remarkable that the value of e c - as estimated 
through the behaviour of 5k or A or by the position of the maximum of k - remains the same 
as in the constant density case, and corresponds to a temperature T ~ 0.15e/fcs , i.e. occurs 
at a temperature of physical relevance for the thermodynamics of this system. Moreover, 
the value of e c is in excellent agreement with the value of the crossover energy estimated via 
the nonequilibrium dynamics fll . 



Let us now briefly comment about the problem of the relevance of quantum effects for 
the results presented in this Section. As already stated in the Introduction, our analysis is 
a completely classical one so that the significance of our results depends on the reliability 
of the classical approximation of the dynamics in the range of temperatures investigated. 
The strength of the quantum behavior of the interaction is ruled by a coupling parameter 
g = hu>Lj/e that is the ratio between the typical vibrational quantum energy Hulj calculated 
in the harmonic approximation and the binding energy e. Values g <C 1 denote that quantum 
effects may be neglected in a wide range of temperatures. However, as observed in Ref. |2B 



the failure of the classical approximation becomes quite evident for all rare gas solids at 
measurable temperatures. The detailed evaluation of the influence of quantum effects for 
the determination of the threshold energy e c is out of the purposes of this paper. We present 
here results for solid Xenon where, due to the small coupling parameter (g = 0.106), the 
relevance of quantum effects in our discussion, if any, is certainly not decisive. 



V. CONCLUSIONS 

We have presented here a detailed description of how a differential geometric approach 
to Hamiltonian chaos can be used for determining the existence of the SST in a realistic 
system. We have discussed the connection between the SST and the geometrical observables 
which have been investigated in this paper. We have shown how standard methods of 
classical statistical mechanics can be used for our computational purposes. We have finally 
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applied the framework here developed to a Lennard Jones crystal modeling solid Xenon. 
The crossover energy region of this system has been clearly detected and the value of the 
energy threshold has turned out to be in excellent agreement with the numerical results 
presented in a recent paper [[□]] . The possibility of setting up an experiment of this system 
for determining the physical consequences of the transition from weak and strong chaos 
regime is under investigation. 
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FIGURES 

FIG. 1. Average Ricci curvature k vs the energy per degree of freedom at constant density 
(p = 1.0 in natural units). Open circles: data of computations where the interaction potential 
(26) has been expanded up to the fourth order in the Taylor series; solid circles: data obtained 
with the "full" interaction potential (26). The two dotted lines represent the estimated low- and 
high-energy behaviours of k helping the identification of e c ; the line interpolating the high-energy 
data has been obtained by fitting a law k ~ a + /?e 1//2 , with fitting parameters a and /3, to the 
high-energy data. 

FIG. 2. Normalized mean-square fluctuation of the Ricci curvature 5^/k vs the energy per 
degree of freedom at constant density. Symbols as in Fig. 1. The dotted line is the linear 
behaviour. 

FIG. 3. Theoretical estimate of the Lyapunov exponent A according to Eq. ( |l6| ) versus the 
energy per degree of freedom at constant density. Symbols as in Fig. 1. The dotted line is the 
power law e 2 . 

FIG. 4. Average Ricci curvature k vs the energy per degree of freedom. Here, in each simulation 
we have used values of the density according to the empirical equation of state of solid Xenon. 

FIG. 5. The same as Fig. 4 for the normalized means-square fluctuation of the Ricci curvature 
5k/ k. The dotted line is the linear behaviour. 

FIG. 6. The same as Fig. 4 for the estimate of the Lyapunov exponent A. The dotted line is 
the law e 2 . 
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